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We present an E(, Grand Unified model with a realistic pattern of fermion masses. All standard 
model fermions are unified in three fundamental 27-plets (i.e. supersymmetry is not invoked), which 
involve in addition right handed neutrinos and three families of vector like heavy quarks and leptons. 
The lightest of those can lie in the low TeV range, being accessible to future collider experiments. 
Due to the high symmetry, the masses and mixings of all fermions are closely related. The new 
heavy fermions play a crucial role for the quark and lepton mass matrices and the bilarge neutrino 
oscillations. In all channels generation mixing and CV violation arise from a single antisymmetric 
matrix. The E(, breaking proceeds via an intermediate energy region with 5*17(3) l x SU{3)rX SU (3)c 
gauge symmetry and a discrete left-right symmetry. This breaking pattern leads in a straightforward 
way to the unification of the three gauge coupfing constants at high scales, providing for a long proton 
lifetime. The model also provides for the unification of the top, bottom and tau Yukawa couplings 
and for new interesting relations in fiavor and generation space. 

PACS numbers: 12.10.Dm, 12.15.Ff, 14.60.Pq, ll.lO.Hi 



The exceptional group Eq P|, 0| is the preferred 
group for Grand Unification. All Standard Model (SM) 
fermions arc in the lowest 27 representation. Its maximal 
subgroup SU{3) x SU (3) x SU{3) can be viewed as an ex- 
tension of the Weinberg-Salam group SU{2)l x U{1)y x 
SU{3)c ^ 5C/(3)l X SU{3)r x SU{3)c = Gasa- The 
fermions can be described by singlet and triplet represen- 
tations of the SU{3) groups only. Using for all fermion 
fields left handed two component Weyl spinor fields, the 
uantum number assignments are (for each generation) 



s 



Quarks 
Leptons 
Antiquarks 



Ql{x) - (3,1,3) , 
L{x) = (3,3,1) , 
Qr{x) = (1,3,3) . 



(■1) 



The 78 generators of Eq consist of the three SU (3) ad- 
joint octet generators and the generators i^(3, 3,3) and 
F(3,3,3) of coset Eq/Gsss- 

The beautiful cyclic symmetry of Eq is apparent from 
(pll and from the fact that F takes a quark field into a lep- 
ton field, a lepton field into an antiquark field and an an- 
tiquark field into a quark field. An additional argument 
favoring Eq is its appearance through compactification of 
the ten dimensional Eg x Eg heterotic superstring theory 
on a Calabi-Yau manifold. The compactification process 
can lead either to four dimensional Eq gauge symmetry 
(which is anomaly free and left-right symmetric) or to 
some of Eq^s maximal subgroups Phenomenology of 
Eq gut attracted attention earlier 0, 0, and its 
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active studies has been continued until recently (6|] . Phe- 
nomenology and properties of G333 triunification models 
are also interesting J]. 

According to (Q) one has besides the SM fermions: ad- 
ditional quark and antiquark fields with the same charges 
as the corresponding down quarks, two SU{2)l doublet 
leptons (containing additional 'active' neutrinos), and 
two SM singlets - 'right handed' neutrinos for each gen- 
eration. 

S'O(IO) and Eq Grand Unified Theories in old 
times usually predicted small neutrino mixings since 
in straightforward applications the large symmetry ob- 
tained from these groups connect the neutrino mixings 
with the small mixings observed in the quark sector. Af- 
ter the observation of large mixings in neutrino oscil- 
lations one had to return to the smaller SU{5) group 
(the minimal version of it does not involve right handed 
neutrinos) or needed several Higgses of the same repre- 
sentation, or special composite operators and fine tuning 
procedures. In this paper we will show, however, that the 
consequent use of the fermion and scalar particle inter- 
actions and spectra of Eq allows to construct a realistic 
GUT model. 

We consider at first the Yukawa sector of Eq with its 
symmetric and antisymmetric matrices in flavor and gen- 
eration space. After defining the model, we can calculate 
from it the mass spectrum of ordinary and new fermions 
and their mixings in terms of a few parameters only. An 
interesting feature is that the mass matrices of quarks 
and leptons are strongly influenced by the flavor mixing 
of the SM particles with heavy fermions as was suggested 
by Bjorken, Pakvasa and Tuan (sl . Earlier suggestions for 
the mixing of the SM particles with new heavy fermions 
can be found in [|| . Our work is done in the spirit of ref. 
0. As in this reference, our scenario favors a relatively 
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light mass scale for some of the new particles [10-plets of 
50(10)]. The lightest can lie in the low TeV region or 
even below. A major difference to 8] is the full use of 
the discrete left-right symmetry of Eq, valid at the inter- 
mediate symmetry G333. It is broken solely by the Ma- 
jorana property of very heavy neutral leptons (the right 
handed heavy neutrinos). The use of an antisymmetric 
Higgs representation proposed many years ago plays 
a decisive role. The corresponding antisymmetric matrix 
determines the generation mixing and the CV violation 
in all heavy and light channels (in the basis in which the 
up quark mass matrix is diagonal). The inclusion of all 
neutral leptons of Eq allows to connect the mass matrix 
of the heavy neutrinos with the diagonal up quark mass 
matrix and the antisymmetric generation matrix. It leads 
to bimaximal mixings of the light neutrinos which then 
changes to a bilarge mixing pattern at the weak scale 
by renormalization effects. All mass ratios and mixing 
angles of light and heavy fermions are simply related to 
each other. 

We then study the gauge coupling and top-bottom-tau 
unification in Eq. It is achieved by an unbroken G333 
subgroup as an intermediate symmetry. The discrete 
left-right 2?l_r symmetry, which is unbroken at these in- 
termediate energies, plays also an important role. The 
breaking scale Mj of the intermediate symmetry is not a 
free parameter, but uniquely fixed by the standard model 
couplings: Mj ~ 1.5 • lO^'^ GeV. Mj also determines the 
scale of light and heavy neutrinos in agreement with ex- 
periment. The unification of the couplings occurs above 
10^^ GeV, in our specific model at 2 • 10"'^^ GeV and thus 
suppresses proton decay. The renormalization of mass 
ratios, various Yukawa matrices and the scaling of the 
neutrino mass matrix are studied in detail. 

Our model is non supersymmetric as the one in 
The hierarchy problem persists but it is hoped that its 
eventual solution would not change the basic features of 
our approach. 



PARTICLE ASSIGNMENTS IN Ea AND THE 
YUKAWA SECTOR 



The charge conjugation operator interchanges left with 
right handed indices: 



C{QL)f = [QRYa , 



CLl C-^ ^L'l 



iQL)f ■ 



(1.3) 



It leaves the commutation relations for the Eq genera- 
tors unchanged and is often called Vlr parity. The new 
lepton fields L\, L|, ig are identical with their own an- 
tiparticle fields. Nevertheless, if two of these fields, say 
L\ and are connected by a single mass term a four 
component Dirac field can be formed. The two fields 
then behave like a (vector like) particle- antiparticle pair. 
The parity and CV operations change the left handed two 
component fields into right handed ones: 

7'(Ql)?(<,x)7'-i =a2(gK)r(t,-x) , 

PL^(t,x)7'-i = (T2Lf (i,-x) , 
V{QR):,{t,x)V-' = ^2(QL)r(^,-x) , 
CV{QL)f{t,x)CV-^ = a2(QL)r(i,-x) , 

CVLl{t,x)CV-^ = a2Ll*{t,-x) , 
CViQR)l{t,x)CV-' = <J2{QRy:{t,~x) . (1.4) 

By including the generation quantum number a{— 
1,2,3) all basic fermions are now classified by the left 
handed Weyl fields '5" with the Eq flavor index r run- 
ning from 1 to 27. 

The product of two 27's of Eq decomposes into a 
symmetric 27, an antisymmetric 351^ and a symmet- 
ric 351s representation 



27 X 27 = 27 + 351^ + 351 



s ■ 



(1.5) 



Consequently, the Yukawa interactions in the Eq La- 
grangian contein in general the three Higgs fields 

H = H{27) , Ha = H{351a) , Hs = H{351s) . (1.6) 

Each of the three Higgs fields couple to the fermions to- 
gether with a 3 X 3 matrix G acting on the generation 
space 



Let us first consider the lowest particle generation 



L\ E- 



E+ Ll 



{Qr)1 = (wa, da, 



(Ll) 



where i,k,a = 1,2,3; a is a color index. In this de- 
scription SU (3) L acts vertically and SU{3)r horizontally. 
The charges are obtained from the operator 

Q = il3 + \y)l + (/3 + \y)r , (1.2) 

with /g, Y defined as usual. Before symmetry breaking 
equivalent forms of (jl.ip can be obtained by applying left 
and right U-spin rotations. 



= [G(27)]„^ , A^p = [G{351a)U . 

So^P ^ [G(351s)]a/3 . (1.7) 



The Eq invariant Yukawa interaction reads 

Cy = ((*J?)^ia2«'f) [Go,f3Hrs + A^piHA)rs + 
Sc,l3{Hs)rs]+h.C. 



(1.8) 



G and S are symmetric matrices in generation space, 
while A is an antisymmetric matrix. Cy is invariant with 
respect to C, the right^left operation. In case of real vac- 
uum expectation values (VEVs) of the Higgs fields, the 
part of Cy obtained from the real part of these matrices 
is formally even under the CV and V operation, while the 
term arising from their imaginary parts is formally odd 
under CV and V. 
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The decomposition of the Higgs fields with respect to 
the G333 subgroup reads 

= (3, 3,1) + (1,3, 3) + (3,1,3) , (1.9) 

= (3, 3,1) + (1,3, 3) + (3, 1,3) + 
(3, 6, 1) + (1, 3, 6) + (6, 1, 3) + (6, 3, 1) + (1, 6, 3) + 

(3, 1, 6) + (3, 8, 3) + (3, 3, 8) + (8, 3, 3) , (1.10) 

Hs = (3, 3,1) + (1,3, 3) + (3, 1,3) + 
(6, 6,1) + (1,6, 6) + (6, 1,6) + 
(3, 8, 3) + (3, 3, 8) + (8, 3, 3) . (I.ll) 

The color singlet parts who's neutral members can de- 
velop VEVs are 

i/^ (3,3,1) , 
ff^^ (3, 3,1) + (3, 6,1) + (6, 3,1) , 

ils (3, 3, 1) + (6, 6, 1) . (1.12) 

We note that the parts containing a sextet or antisextet 
representation can only couple to leptons. 

II. THE MODEL 

The vacuum expectation values of the three Higgs 
fields determine the particle spectrum. To be in accord 
with the SM the masses of the new particles of Eq have 
to get heavy (at least of order TeV). Thus, Higgs com- 
ponents which are SU{2)l singlets can have large VEVs. 
The members of SU{2)l doublets, on the other hand, 
should be of the order of the weak scale, while the VEVs 
of SU{2)l triplets are expected to vanish. 

In order to define our model to be predictive and to 
have very few unknown parameters, we need some spe- 
cific assumptions concerning the three Higgs fields, about 
the generation matrices G, A, S and the symmetry break- 
ing pattern. We do not consider Higgs field components 
which carry color. They are supposed to aquire masses 
of the order of the GUT scale from appropriate Higgs 
potentials. 

We allow VEVs for all color singlet and neutral com- 
ponents of H 

{Hi) - e\ , {HI) = el , for i, fc = 2,3 . (H.l) 

However, by a biunitary left and right U-sjim transfor- 
mation in flavor space (i.e. on the SU{3)l and SU{3)u 
indices 2 and 3) we can choose a proper basis for which 

el = el^ 0. (H.2) 

Our first assumption concerns the VEVs of Ha and Hs- 
{Ha) can mix the standard model particles with the new 
heavy D and L states [the 10-plet of 5*0(10)]: d ^ D, 
e E ., V = L"^ L\. This is achieved by components of 



Ha which involve left and right W-spin 1/2 indices. For 
the (3, 3, 1) sector of Ha we take, therefore, 

{iHA)l)^fl, i,k^2,3. (H.3) 

For the (3, 6, 1) sector of Ha one has correspondingly 

{Hf-'^) ^ r^'^'K *,A: = 2,3, (H.4) 
and for the sector (6,3,1) 

{HA{i,k}i) = /{i,fc}i , i, fc = 2, 3 . (H.5) 

In our numerical treatment we will restrict the VEVs in 
(dU, (P3|l to those with i = 3, A: = 2, 3 and i = 2, 
fc = 3 which should be the dominant ones. With respect 
to ht-spin, /•^f i^'^l is the analogue of /I- While /I mixes 
d with D, /•^f i''^! mixes e~ with E~ and v with L|. 

The VEVs of the symmetric Higgs field Hs can provide 
large Majorana masses for the heavy leptons L2 and L^. 
They arise from the Hs{6, 6, 1) sector. Here we have to 
take the SU{2) l singlets and left and right handed U-spin 
triplets 

All other components of {Ha) and {Hs) are taken to be 
zero or negligible in our calculations. 

Of particular interest is the question of the breaking 
of the left-right symmetry of Eg. F^^-^^ as obtained from 
{Hs) breaks this symmetry strongly. It could be the dom- 
inant manifestation of X'Lfl symmetry breaking. {H) and 
{Ha) on the other hand need not break this symmetry 
significantly. A strict left right symmetry in this sector 
would imply the relations 

fk = -ft : r^'-'^ = , ^,k = 2,3. {11.7) 

The signs follow by taking the Ha part of the Yukawa 
interaction to be even under Vlpi when Ha is replaced 
by {Ha)- As a consequence of (III.7|) the f's are of the 
order of the weak scale even though some are standard 
model singlets and thus only protected by the discrete 
T^LR symmetry itself. If this is indeed the case, it implies 
new particles in the few TeV region as we will see. 

The next assumption concerns the generation matrices 
G, A and S. The symmetric matrix Gap can be diago- 
nalized by an orthogonal transformation, which leaves 
the symmetry properties of Aap and Sap unchanged. By 
choosing this basis, the up quark mass matrix is diago- 
nal because, according to the above assumed properties 
of {Ha) and {Hs), only {H) contributes to it 

i'^u)a0 = Gapel = gaSafsel . (II. 8) 

As a consequence, the quark mixing angles and the CV 
violating phase must entirely come from the inclusion of 
the Higgs Ha with its antisymmetric generation matrix 
Aap as proposed in ref. 0. Thus, Aafs has to contain 
imaginary parts which can not be rotated away using 
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quark phase redefinitions. This leads us to assume that 
the matrix A is - in our phase convention - purely imag- 
inary, i.e. a hermitian matrix. The normalized matrix 
contains then only two parameters, in fact only one when 
utilizing a discrete generation exchange symmetry for A 
as shown later. 

We suggest, that the generation matrices G, A and 
S are not independent of each other. In particular, the 
coupling matrix S for the heaviest leptons should have 
an intimate relation with the generation matrices of the 
charged fermions S may then be expanded in terms 
of G and A. Speculatively we assume: The generation 
mixing matrix S* is a combination of the bilinear product 
G^ and the commutator [G, A]. The generation mixing in 
this sector is then due to the same matrix A which causes 
the mixing of the charged fermions. As it turns out this 
structure for S is crucial for bilarge neutrino mixings. In 
fact, it leads to bimaximal mixing which is then changed 
to bilarge mixing by renormalization group effects. 

The last assumption concerns the breaking pattern of 
Eq , which is presumably the origin of the breakings seen 
in the Yukawa sector. We suppose the following symme- 
try breaking chain: 



and lepton masses at the scale /j, = Mz 



ruu = (1.8 ±0.4) MeV , 
m, = (62 ± 12) MeV , 
nib = (2.89 ±0.03) GeV , 
me = 0.487 MeV , 



md = (3.3 ±0.7) MeV , 
mc = (0.64 ± 0.04) GeV , 
mt = (173 ± 5) GeV , 
= 102.8 MeV , 



nir = 1.747 GeV , 



(11.10) 



as obtained from the analysis of experimental data. The 
general hierarchical structure of the SM masses and of 
the CKM matrix elements will be used in the following. 
Some of the masses, in particular mt, rub and rrir, are 
taken as input parameters. 



III. THE QUARK MASS MATRIX 

Because of the hierarchical structure of the quark 
masses and mixing angles it is convenient to express them 
in terms of powers of a small dimensionless parameter. 
We introduce the parameter a with the value 



a = 0.058 



(IILl) 



for which 



El 



Mgv 



Mi 



' SU{3)l X SU{3)r X SU{3)c x Vlr 

SU{2)l X U{1)y X SU{?.)c ■ (11.9) 



Here Mqut is the GUT scale and I?l_r denotes the dis- 
crete left<-!-right symmetry operation. As we will show 
below, the breaking chain ljll.9|) leads in a straightforward 
way to the unification of the gauge coupling constants. 
The first breaking step to the intermediate symmetry can 
be caused by a scalar 650-plet which contains two G333 
singlets. One of them 5+ is even under Vlr (5+ 5+), 
while the second one S- is odd [S- It then 

follows from the symmetries at Mj and above Mj that 
S+ has the non zero VEV and = 0. This insures 

that in the {Mj, Mgut) interval L ^ R symmetry is 
precise and the equality of the coupling constants gL and 
gij is protected also at the quantum level. 

We take two Higgs SU (2) l doublets of H, namely ff^'^ 
and H2'^ , to be relatively light. The remaining Higgs 
masses of the color neutral components of H and Ha are 
taken to be of order Mj or higher. The only exception is 
the SU{2)l doublet {Ha)2 which can be much lighter 
than Ml because of the left right symmetry (/| « 0) in 
the Ha sector mentioned above. But it must be heavier 
than « 500 TeV not to induce flavor changing processes 
above presently known limits. All Higgses not mentioned 
are assumed to have masses at the order of the GUT 
scale. 

Before starting our investigation, let us state the quark 



\Vcb\ ^ and 



mt 



cr^ (IIL2) 



holds within experimental uncertainties. One also has 

\ ■ to our assumption ljll.l|) , ljll.2l) the up-quark 
mass matrix is 



Diag f 2hL, ± 1) mt . 

At the scale ~ Mz we can write 

mu Diag (cr^ ± cr^, l) 173 GeV . 



(IH.S) 



(in.4) 



The signs of the mass parameters are in general of no 
relevance because of the freedom to change phases. But 
since we keep G and A to be hermitian matrices, the 
Jarlskog determinant obtained from the commutator of 
mass matrices depends on the sign chosen in (|111.3|) giving 
two solutions for the area of the unitarity triangle. 

Because of the existence of the D-quarks, the down 
quark (big) mass matrix is a 6 x 6 matrix, which contains 
the antisymmetric generation matrix A 



Md,D = 



d 
D 



\G 



D 



elG^ 



fiA 



(111.5) 



Here 63, /| and /| are mass scales of order of the weak 
scale, while at least 63 should describe a heavy mass scale. 
In accord with our model assumptions we have /| , /| ^ 
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e| which allows to integrate out the _D,_D states and to quarks and leptons. Below Mj, one gets 



write down the see-saw formula 



niD ^ elG + f^A 



fifi 



(AG-^A) 



(III.' 



The D-quark mass matrix is simply proportional to the 
up quark mass matrix: 



D 



(III.7) 



Although eq. (|III.5I) should only be valid at the unifi- 
cation scale and has to be carefully scaled down for a 
determination of m^i at /i ~ Mz, we will use eq. Ijlll.6|) 
at Mz for a first orientation. 
The first entry in 



l|III.6|l is responsible for the mass of 
the bottom quark, while the second term must provide for 
the small mixing angles and the large CP violating phase. 
The third term gives a correction to the symmetric part 
of the mass matrix which is important for the strange 
quark mass. We expect, therefore. 



' nib 



23 



\Vcb\mb 



/||^i2| ~ m,\Vus\ , and /||Ai3| ~ IKolm, (III.8) 

From ljIII.2|l one then gets for /| and the generation ma- 
trix A 





( »• 


icT , 


—id 




-icr , 


, 






^ iff , 


"71 ' 






(III.9) 

We introduced a scaling factor \/2Xa (as discussed in sec- 
tion such that /I = and the scale dependent 
matrix A is normalized according to Tr(A^) ~ 2A^. We 
remark that the antisymmetric matrix A taken here is 
also antisymmetric with respect to the (discrete) inter- 
change of the second generation with the third one. We 
know of course, that the matrix A can have its strictly 
antisymmetric form only above M/, the breaking point 
of the left-right symmetry. Thus, in our renormalization 
group treatment we take the matrix A as given in eq. 
(|III.9|I to be strictly valid at ^ = Mgtjt , even though we 
anticipated its form at a low scale. As we will discuss in 
the appendix, by going down from fi — A/gut to M/, the 
matrix A 'splits' into a matrix A'^ for the quarks and a 
matrix A^ for the leptons. By going further down to Mz, 
A'^ as well as A^ each splits into three matrices relevant 
for the sectors indicated by the superscripts: 



AQ 

A^ - 



+ 



Dd 



and 



A 



e e ' A^ A^J e+ 



(III.IO) 



These matrices are no more strictly antisymmetric. Ob- 
viously, also the matrix G splits into more matrices. Be- 
tween Mqut and Mi we have G (Gq, Gl) for the 



Q — ^ 1 ^ 7 ^ 7 ^ 

Gl (g"^"' , G^^ , G""^ 



DD 



and 



(III.ll) 



We calculated these matrices at /i = Mz in a model 
specified in section IVlII which has a unification scale of 
2 • lO^'^ GeV. The matrices A are exhibited in (|A.45p - 



The matrix 



becomes e^G''''. 



It is 



obtained from mu = ei[G"" at Mz [eq. 1111. 4() ] scaled 

up to the GUT scale, where G"" = G'^'^ holds and then 
scaled down to Mz- In our approximation, it is still 
a diagonal matrix. In the same way, one obtains the 
matrix replacing G^^ in the DD channel of eq. l|III.6p . 
It is denoted by {G^^) -\ Th e matri ces G^^ and G'^^ 
at /i = Mz are given in HA.12p . HA.33|) . 

With these changes the mass matrix for the down 
quarks becomes now 



moiMz) 



elG 



dd 



flA'^ 



(G^^)-M^'^ 



b ' 
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^/2 



(III.12) 



rUb will slightly differ from the mass of the bottom quark 
because of the mixing occuring in mo ■ 

After having found the renormalization group effects 
on the matrices G and A, the only parameter for cal- 
culating the d-quark masses and the CKM matrix is 
fifi/{el{G^^)33). We use this parameter for a fit of 
the Cabibbo angle \Vus\- Because of our expectation of 
an approximate left right symmetry [see (|II.7p ] we look 
for a negative value of this parameter and find 



fif2 



^(G^^) 



-4.75 • 10-^ GeV 



(III.13) 



33 



Upon diagonalization of the down quark mass matrix 
with the negative sign taken in ljIII.4p . one ob- 



md{Mz) ^ 2.66 MeV , ms{Mz) 
nibiMz) ~ 2.89 GeV , 
|Ks|=^ 0.217, [1461^ 0.045, \Vub\ 



49.7 MeV , 



0.0034 ,(111.14) 

and for the angles of the unitarity triangle 

a -84°, ^-20°, 7 -76°. (III.15) 

To obtain the correct value for mb{Mz) we took for the 
(3,3) element of e^G^Vn^ = 2.859 GeV. A similar good 
fit is obtained if in ljIII.4|l the positive sign is chosen. The 
number given in HIII.13|) then changes to -3.26-10-5 GeV 
and the angles of the unitarity triangle become 



a ' 



95°, 13-21° 



7 ~ 64° 



(III.16) 
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In the following we will use the negative sign in (|III.4|) . 

The results H111.14ll -( |111.16|) are in good agreement with 
present experimental data. The mass of the strange 
quark is a bit low but still within the bounds of (jII.10|) . 
We also see, that Weinberg's suggestion 

(III.17) 

is valid. It follows from the smallness of the (1, 1) entry 
62(7'* in (|III.12p due to the small first generation up quark 
mass. We further note, that the term in mu, which arises 
from the mixing with the heavy D-quarks, reduced the 
angle 7 from the originally obtained value ~ 90° 5] to a 
lower value. 

Besides (IIIL13|I there is no restriction on the value of 63 
except that fi/e^ has to be sufhciently small to justify 
the see-saw formula and thereby the near unitarity of 
the CKM mixing matrix. However, as mentioned in sect. 
2, the VEVs (H) and (Ha) may approximately respect 
the left-right symmetry of Eq and of the intermediate 
symmetry in contrast to the large VEV of Hs ■ This idea 
is supported by the small value found for /| in l|III.9(l . 
It would be zero for a strict left right symmetry in this 
channel and is indeed small (/| ~ 0.093 GeV) compared 
to the weak interaction scale. One can then expect, that 
the product — /f/l is not of order AIzMj but not much 
higher than (Af^)^. This gives us a rough estimate for 
63 and thus for the masses of the D quarks. 

e3(G )33 ~ -Q^^Jz or ^ -^^M, . 

(III.18) 



From these relations, which are of course sensitive to 
the value taken for the weak scale input, we expect 
4iG"^)33 to be of order (10^ - 10^) GeV. Taking 
63 (6*^^)33 = 4 • 10^ GeV as an example (and scaling 
effects into account), one obtains 

Mdi =i 557 GeV , Md2 ^ 129 TeV , 

Md3 ~ 4 • 10^ TeV . (III. 19) 

A more detailed discussion of the heavy fermions and 
their masses will be presented in section IVII and in the 
appendix Al. 



IV. THE CHARGED LEPTON MASS MATRIX 

The charged lepton mass matrix has the same structure 
as the down quark mass matrix. By going from quarks to 
leptons Eq Clebsch- Gordon coefficients have to be taken 
into account. Quarks and leptons couple to H(3, 3, 1) 
according to the combination 

q^q' + ler,."e''''^"LiLt:, . (IV.l) 

The (3, 3, 1) sector of the Higgs field Ha couples only to 
quarks, the sectors (3, 6, 1) and (6, 3, 1) only to leptons. 
Thus, the relevant 6x6 matrix at the GUT scale is 



Me.E = 



E- 



-/{l,3}3^ , 



E+ 

/3{l,3}^ 



(IV.2) 



Using the same arguments as for the down quark mass 
matrix the /'s in the diagonal elements are small com- 
pared to the main terms. After integrating out the E- 
type states, the mass matrix for the charged leptons of 
the SM is generated and has at /i ~ Mz the form 

ruE =^ -elG^-'^^ - (/2{i3} _ f^,,^,)A^--^ 

The first term is constructed like e^G'^'^, but for leptons 
and given in the appendix A2. The contribution of VEVs 
in the second term should be as small as the correspond- 
ing term /| in the quark mass matrix. Diagonalizing 



(jIV.3p one gets with 

/'^"^ - /{13}2 ^ 0.042 GeV and 
^^'iV^lll''"^ - 12.6 • 10-5 GeV , (IV.4) 
the charged lepton masses 

TOe = 0.488 MeV , m^, = 102.8 MeV , = 1.748 GeV. 

(IV.5) 

For obtaining the correct value of the tau lepton mass we 
took m° [the (3,3) element of elG"'"^] to be 1.6 89 GeV. 
The contributions from the first term in ljIV.3|l for the 
light generations are proportional to ct'* and a'^, respec- 
tively, and thus negligeably small. The muon mass re- 
ceives its essential contribution from the third term in 
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(|IV.3|I . i.e. from the mixing with the heavy leptons. The 
contribution from the second and third terms to the elec- 
tron mass are comparable. There is some CT'-violation 
due to the second term in (|IV.3|) . The corresponding 
unitarity triangle, for charged leptons, has the angles: 
a ~ 43° , /3 ~ 62° , 7 ~ 75° . After diagonahzation of 
the charged lepton matrix, this CV violation will affect 
the neutrino mixings. The charged lepton mixing an- 
gles turn out to be small |V;^| ~ 0.034, \Ver\ ^ 0.003, 
\V^t\ — 0.068. Therefore, the large neutrino mixings 
are not due to the mixings in the charged lepton sector 
but should come from the neutral lepton sector where 
large Majorana masses appear. In the next section it 
will be shown that this is indeed the case. Comparing 
now l|III.13(l with (|IV.4|I and taking into account ljA.18|l 



we get ^^'^'^La — (1-6)^. Considering the analogy 

of /| with /{i,3}3 and /| with /^li'^} ^j^jg appears to 
be a reasonable value. In sections IVII IVIII we will use 
l/{i,3}3l ^ l/ll - l/ll, i-e. appropriate left- 

right symmetry in these channels. 



V. THE NEUTRAL LEPTON MASS MATRIX 

The fundamental fcrmion representation of Eq contains 
five neutral two-component fields. Thus, for three gen- 
erations, the mass matrix for these neutral leptons is a 
15 X 15 matrix. According to the assumption stated in 
section m it is given by 





Ll 


Ll 


L\ 

_/3{l:3}^ 


Ll 





-e\G 







-e\G 







~/{l,3}3^ 

elG+hlA 













e\G 


/3{1:3}^T 


-/{l,3}3^^ 







elG+hlA 








e\G 


elG+hlA^ 


/ 



where 

hi = f^'-'^ + /{1,3}2 , hi = ^ f^^^^^^ (V 2) 

and ig stands for the standard light neutrino fields. All 
ingredients in this matrix arising from the Higgs fields H 
and Ha are defined in the previous sections. We notice, 
however, that in the L\Ll block the contribution of the fs 
is additive. The new elements are the ones containing the 
symmetric generation matix S. They give rise to genuine 
Majorana mass terms and are of particular significance 
in the diagonahzation process. The strength of the Hs 
Higgs contribution to is governed by the constants 
^{2,2} ^j^^ ^{3,3} carrying right handed U-spin quantum 
numbers. F^^''^'^ essentially fixes the Majorana mass for 
the heavy Lg leptons which are expected to be of the or- 
der of the SU (3)l X SU (3)i^ breaking scale. The constant 
^{2:2} similar strengths breaks the left-right symmetry 
and thus is responsible for the dominant breaking of this 
symmetry. 

We can reduce the matrix Ml to a 9 x 9 matrix by 
knowing that e| is much larger than the other elements 
in the same row and column, in particular, if /-^{^^■^l and 
y3{i,2} g^^.g indeed of the order of the weak scale. This 
allows to integrate out the L\, Ll states. With the ab- 
breviations 

/ = /3{i,3}eJ/e3 , / = /{i,3}3e}/eL (V.3) 



one finds 

Ll / -e\G fA \ 
Ml' = Ll ~e\G F^^^^^S fA , (^.4) 
Ll\fA^ JA^ Fi^^^^sJ 

and 

MlilI - 4G ■ (V.5) 

We neglected in (jV.4|) a correction to the (3-3) block, 
namely —-^e2G. It is small compared to the large eigen- 
values of F^^'^^S. 

For values of /3{i^3} ^j-^g qj-^^j- of the weak scale 
and F«2.2)}^^{(3,3)} jjgj^j. ^ iQi3 Qgv, 

we can again 

apply the see-saw mechanism and finally arrive at the 
3x3 Majorana matrix for the light neutrinos 

^ -^GS-^G - j^A-S^^A , (V.6) 

and for the mass matrices of the heavy Majorana neutri- 
nos 

Mli ^ i^{2,2}_5. ^ ^^^^ _ ^{3,3}^. (Y.7) 

Only the first term in (|V.6p need to be considered, since 
the remaining one can safely be neglected. Therefore, the 
neutrino mass matrix (|V.6p is mainly due to the decou- 
pling oi Ll = V states. It scales with the masses of these 
heavy lepton states. 
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We expect [see sect. 2] S to be related to the two 
other generation matrices G and A. The main term for 
S should be G^, which leads to a diagonal non degenerate 
mass matrix [see (|V.6|l ]. We then add a term linear in A, 
the commutator [G, A], with a tiny coefficient. It implies 
that also in this sector generation mixing is solely due to 
the antisymmetric matrix A. We take for S, divided by 
the overall coupling strength As to the Higgs field Hs, 
the real and bilinear construct 



S/Xs = G^ +ixa^[G, A] 



(V.8) 



whith the single parameter x. The G^ term with its dom- 
inant element ~ 1 for the 3"'^ generation serves for gen- 
eration hierarchy and for the normalization of S/ Xs (for 



which the term in (|V.8p can be neglected). With no 
renormalization effects included, the matrix S, as defined 
in IjV.Sp reads 



S/Xs 



( 




-a^x 
















4 


^x 


V 


—ax , 





71 " 



(V.9) 



In each element of (|V.9|I only the leading powers of a are 
shown. 



By inverting the matrix S defined in ljV.8|) and using 
(jV.6p . one finds for rrii, 



( -1 



(1 



^ax)x 



- 1 , 



-(1 

{x- 



^ax)x 



(l-2£!;2-|- Vaster) As-F{2, 2} 



(V.IO) 



For the simplicity of representation ljV.10|) contains only 
the zeroth and first powers in a. Taking the full expres- 
sion makes numerically little difference. The interesting 
feature of is the fact that it produces for any value 
of a; ^ 1.5 automatically an almost perfect bimaximal 
neutrino mixing pattern (!) with a normal (not inverted) 
neutrino spectrum. By changing a;, solely the ratio of 
mass square differences 



R 



(V.ll) 



changes {rrii denotes the three ly eigenstates mass ordered 
according to mi < m2 < m^). The experimentally ob- 
served ratio R « 0.03 is obtained for x ~ 3.5. 

However, for a proper calculation of the neutrino mass 
matrix at fi = Mi and fi = Mz, renormalization effects 
have to be taken into account. This is particularly neces- 
sary because of the large generation splitting of the heavy 
neutrino states L2 — v caused by the G^ term in the ma- 
trix S. We have to integrate out these states in steps 
and to redefine m,„ in each step. We start by using HV.8() 
at the scale Mi with G = Gl and A = Al and proceed 
according to the rules given in appendix A4. It turns out 
that renormalization effects strongly infiuence the neu- 
trino mass matrix and thus also the mixing pattern. The 
bimaximal mixing is changed to a bilarge mixing. The 
calculation is again performed for the gauge and Yukawa 
unification at 2 • 10^^ GeV described in section IVlII As 
in the examples given in [Tsl l we find that the renormal- 
ization coefficients strongly reduce the mixing angle 612 
observed in solar neutrino experiments while the angle 



023 observed in atmospheric neutrino experiments is less 
affected. The renormalization coefficients also increase 
the value of the ratio R. 

A good description of the known neutrino data is ob- 
tained by changing the value of our parameter x = 3.5 
to 



2.8 



(V.12) 



With this value we obtain R ~ 0.055. Larger values of 
X reduce R. However, this would lead to a too strong 
reduction of the solar neutrino oscillation probability. 

With X ~ 2.8 one obtains for the mass matrix of the 
light neutrinos at /i = Mz 




Mo (V.13) 



with 



Mo 



(A^(M/)ei(M/))2 (102.17 GeV)2 



Me 



(V.14) 



Here A,- is the coupling of the third generation lepton to 
H2 '^ and e\ is the VEV of the Higgs field as used 
before. 

We obtain the mass squared difference observed in at- 
mospheric neutrino experiments Arn^^^ ~ 2.5 • 10-3 eV^ 
when setting M^^c^ ~ 1.6 • 10^* GeV. It is very satisfy- 
ing that this scale is of the same order of magnitude as 
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expected from the value of Mj, the breaking point of the 
left-right symmetry With this value of M^^p^ the 

neutrino mass eigenvalues turn out to be 

mi = 0.0023 eV , ma = 0.0120 eV , mg = 0.0516 eV , 
m| - m? ~ 1.4 • lO-"* (eV)2 , 
ml - ml ~ 2.5 • lO^^ (eV)^ . (V.15) 

To obtain the neutrino mixing matrix, one has to go 
to a basis in which the charged lepton mass matrix is 
diagonal. Diagonalizing ljIV.3|l and denoting by h'a the 
weak eigenstates {a = e, /i, t), we find from (jV.lSp 

l^c = U^^'y^ , (V.16) 



/ -0.2 + 0.861, -0.46-0.08i, 0.026 \ 
U^,~\ -0.29- 0.0151 , 0.044 + 0.551, 0.78 

V 0.37- 0.0171, -0.036- 0.691, 0.62 / 

\ /at 

(V.17) 

Here we took a special phase choice for the neutrino flavor 
eigenstates such that the 3'''^ column has only real and 
positive elements. 

Our results for the three mixing angles relevant in neu- 
trino oscillation experiments as obtained from eq. ljV.17|) 
are 

sin^eiia ~ 6.8-10-'^ , sin^ 6*12 ~ 0.22 , sin^ 6*23 =i 0.61 . 

(V.18) 

These values and the ratio R ~ 0.055 of mass squared 
differences are quite close to the result of Su- 

perKamiokande fjl , SNO ^9] and KamLAND 
the CHOOZ limit |22| and the observations of the disap- 
pearance of solar neutrinos ^231 . 

We can also get from (|V.17|) the neutrino unitarity 
triangle, defined in analogy with the quark unitarity tri- 
angle. It turns out to be: 

a, ~ 74° , (3,-6° , 7, ~ 100° . (V.19) 

The phases of the elements of the first row of are 'Ma- 
jorana phases' relevant for neutrinoless double /3-decay 
experiments. With the convention used in l|V.17|l we get 

^11 ~ 103 , 612 ~ -170° , 5i3 = . (V.20) 

For the the quantity {m,)ee which determines the decay 
rate we find 

(m,)ee = |mi(C/,^i)2 +m2(t/,^^)2 +m3(;7,''3)2| ~ 

9.5 • 10-* cV . (V.21) 

The matrix Lf^ , in particular its deviation from bimax- 
Imal mixing, depends via the renormalization parameters 
to some extent on the way unification is obtained. But 
a bilarge mixing with near maximal mixing in the fJ. — t 
sector will always result from the basic assumptions of 
our model outlined in section Hll 



VI. THE DESERT IS BLOOMING 

In our model the masses of the heavy down quarks D 
and the corresponding leptons L, which form a 10-plets 
of 5*0(10), have a generation splitting similar to the up 
quarks. The absolute values of these masses can not be 
given. However, if {Ha) still respects to some extent the 
left right symmetry of Eq as discussed above, the lightest 
D and L states lie in the TeV region. In section IVIII we 
present a numerical solution of the problem of the gauge 
and Yukawa coupling unification for A/gut — 2-10*^ GeV 
and el{G°^)33 = 3.98 • 10^ GeV. This solution also fixes 
the so far undetermined VEVs of H and Ha- With the 
values quoted there [eq. ljVII.361) ]. we can now directly 
diagonalize the 6x6 matrices (|III.5|) and (|IV.2() which 
determine the mixing of the SM particles with the D 
and L states. This mixing, although important for the 
mass matrices, does not seriously violate the unitarity 
relations for the SM particles. For example, the sum of 
the squares of the second row of the CKM matrix differs 
from one only by 1.2 • 10^^. In the charged lepton sector, 
the corresponding deviation amounts to 3.7 • 10~^. 

We can list now the mass values of the new particles 
by using again 63(6*^^)33 = 4 • 10'' GeV and setting 
A5F{2.2} ^ AsF{3^3} ^ i Q . iQi4 accord with the 

neutrino results: 

Md, 557 GeV , M^, ~ 129 TeV , 

Mzj3 -4-10'' TeV , 
Ml, 355 GeV , Ml^ ~ 103 TeV , 
~ 3.83 • 10^ TeV , 
M{Ll)i ~ 5.8 • 10^ GeV , M(Li)2 ~ 9.6 • 10^ GeV , 

M(L3)3 ~ 1.6 • lO^-* GeV , 
MiLl)i ~ 5.8 • 10^ GeV , MiLl)2 ~ 9.6 • 10^ GeV , 

Af (Li)3 ~ 1.6 • IQi-* GeV . (VI.l) 

In the evaluation we took the most important renormal- 
ization effects into account (see section IVHI and the ap- 
pendix). As we see, the desert is populated between the 
mass scales Mz and AfcuT- The mass ratios for different 
generations of the standard model singlet neutrinos are 
even more drastic than the corresponding ratios for the 
D quarks and the SU{2)l doublet heavy leptons. 

Our specific unification model allows to calculate nu- 
merous properties of the old and new particles in par- 
ticular those related to their decay properties. We will 
present here a few examples only. 

From the 6x6 mass matrix (|III.5ll . for the quarks, one 
can calculate the coupling matrices in generation space 
for the couplings of the light and heavy mass eigenstates 
to the appropriate light Higgs field components 

102^^ d Hi , 1(72^^^ D Hi , 

ia2C^'^ d Hi . (VI.2) 

We find, without using the remaining freedom of chang- 
ing phases, 
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(-0.18 + 1.61) 
(7.4 + 0.28 i) • 
-(9.9 + 0.4i) 



• 10-4 
10-4 , 
10-3 , 



(-0.25 + 7.31) • 10-4 
(3.3 + 1.3 i) • 10-3 , 
-0.043 + 0.0141 , 



-(3.2 + 0.91) • 10-5 
(-0.42 + 1.51) • 10-4 
0.069 + 1 



, 

-(0.3 + 3.41) • 10- 
-0.77-8.91 , 



-4.4- 10" 
, 
0.39 , 





4.11- 10" 




10- 



C 



Dd 



0, (0.33 -3.81)- 10-2 , 0.84-9.71 
-4.8-10-^, 0, 0.43 

0, -4.6-10-6, 



10- 



(VI.3) 



Of course, similar results can be derived for the lepton 
couplings to the Hlggs fields. 

For the weak Interaction process D ^ u Wl one can 
Introduce the matrix V"^ as an extension of the CKM 
matrix 

IZ7^(1-75)V"^D(Wl)+ . (VI.4) 
From (|III.5|) one gets 

/ -3.5 -10-4, -0.04, -1.1- 10-4 \ 
V"-D~ -0.95+111, 0, 1.3-10-3 -10-3. 
\ 0.81-9.41 , 0.42 , j 

(VI.5) 



There are also right handed current Interactions of the 
standard model particles with the heavy SU{2)b, vector 
bosons 

u7^(l+75)V"'^"d(Tyfl,)+ , (VI.6) 

where V"'' Is slightly different but has the same structure 
as the CKM matrix. 

Of particular interest for the decay properties of the 
mass elgenstates of v neutrinos are the Dlrac masses con- 
necting the flavor elgenstates of the light neutrinos (in 
a basis in which the charged lepton matrix is diagonal) 
with the heavy neutrinos. Using (|IV.3|I . e\GL{Mi) from 
(|A.8|I . HA.74|I . HA.75|) and dlagonallzing S{Mi) we obtain 



Vi V2 fs 

vi ( (1.7- 101) - 10-4 , (-1.2 -111)- 10-3 , -0.41 - 0.00991 

m^^^^ ~ (-4.6 + 0.31) - 10-3 , 0.033 + 0.321, 3.8 + 7.71 jCeV. (VI-7) 

vA (3.3 + 3.91) - 10-3 , -0.034 + 0.0411, -93.7 + 83.61 



VII. UNIFICATION OF COUPLINGS 

A. Gauge coupling unification with intermediate 

517(3)1, X SU{Z)r^ X S'[/(3)c symmetry 

As it is known, the SM does not lead to the unifica- 
tion of the gauge coupling constants. In our scenario, 
there are the additional Dlrac fermlons D and L below 
the GUT scale Mgut- However, these do not alter the 
unification picture of the standard model significantly. 
We still need to introduce an intermediate breaking scale 
Mi. 



A large group like Eq with high dimensional repre- 
sentations should first be broken by a step which lowers 
the symmetry considerably. It is natural to break Eq to 
the maximal subgroup SU{i)L x SU{'i)R x SU{i)c- As 
we will see, this has the advantage that the correspond- 
ing intermediate scale is not an arbitrary parameter but 
fixed. The breaking at the GUT scale can be achieved in 
the scalar sector by a Hlggs -ff(650), which contains two 
G333 singlets (1, 1, 1), 5+ and S+ is even under Vlh 
and thus keeps the left-right symmetry, while 5_ is odd. 
We have to take (5_) = and (5+) to be different from 
zero for the breaking. It keeps gL = 5fl [dL = 5s;7(3)t: 
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gR = 5sc/(3)h] for ^ Mj. The reason is, that at the 
intermediate scale M/ the SU{2)l gauge couphng 32 (m) 
and the hypercharge coupling 51 (/z) have to respect the 
SU{?,)l X S'C/(3)i? symmetry. Since the U{1)y hyper- 
charge is a combination of 1^, Hr ^^nd Yr, according to 
(|1.2p . the intermediate symmetry automaticaUy requires 
the matching 



and 5rL(M) = for 



'i)^gL{Mi)=gR[Mi) 

' - > A//. (VII. 1) 



The relation gL{l^) — for M - holds even at the 

quantum level since it is protected by V^ji parity. As a 
consequence, Mj is fixed by the meeting point of 172 and 
gi . From thereon the two curves continue as a single one 
up to Mgut where gL = gn. unifies with gc = gsu(3)c- 
For this to happen the states i?A(6, 3, 1) and Ha{3, 6, 1) 
will play a central role as we will see shortly. 

The details are as follows: Below AIj, the field con- 
tent consists of the fcrmionic generations of the standard 
model together with two light Higgs doublets and the 
three Dirac particles D{D, D), L^{L\,Ll), L-{L\,Ll). 
The two Higgs doublets are hI'^ , H^"^ with (Hl)'^ + 
(H'i)^ = i;2 f ^ (174 GeV)2. vo wiU be smaller than 
V in case an additional Higgs meson with standard model 
quantum numbers has a non zero VEV. 

The corresponding 6-factors for the evolution of the 
couplings are 



21 



-3,-7 



(VH.2) 



as obtained from the standard model fermions and two 
Higgs doublets. The additional 6-factors for the D's and 
L's for each generation are 



(61,62,63) 



D 



4 2 
15'«'3 



(61,62,63) 



2 2 

5'3 



,0 



(VH.3) 

Apart from these additional states, there are more scalar 
doublets i/™ (m = 1,2,3,4), which are involved in the 
construction of the fermion sector. One of them comes 
from _ff(3,3, 1) and the other three from _ffyi(3,3, 1). H 
and Ha also contain two isosinglet fields H" [n = 1, 2) 
carrying the same U{1)y charge as . If some of the 
corresponding states with masses ^{H^) and ^{H^) he 
below Mj, each of them will contribute to the 6-factors 
according to 



(&l,^2,&3 



1 1 



,0 



(^1,^2,^3) 



^4 



1 



:,0,0 



10' 6' y ' ^ - - - \^5' 

(VII.4) 

Four more Higgs components which are SM singlets could 
also be relatively light, but they do not contribute to the 
running of the gauge couplings. Thus, the solution of the 
renormalization group equation (at one loop level) for the 



gauge couplings at M/ reads 

2n l^m^^-iKHf) 2^ TW^MhT) 



■In- 



M, 



K 1 M, 

■ m ■ 



(vn.5) 



Here 6^, Md and 6^, Ml denote masses and 6-factors of 
D and L states respectively. The matching of gi and g2 
at Ml gives 



In^ 

Mz 



-skill- 



\Mz 



Mi 



^(Mz)) + TME„.ln 



>7 A/r,, Mr,„Mr,„ " 1.^11. u; 



Ml 



At the GUT scale we should have Ml, = Md,- Ac- 
cording to sections IIIII to Ml^ — Mjj. should hold 
approximately also at lower scales, since they are deter- 
mined by (H^). Thus, for the determination of Mj we 
can safely neglect the last term in (jVII.Gp . Taking the 
masses /i(iJ™) ~ l^iH"^) ~ M/, also the second term can 
be neglected. With a~[^{Mz) = 59 and a:^^{Mz) = 29.6 
we then obtain for M/, the breaking point of the inter- 
mediate symmetry, Mj ~ 1.3 • 10^^ GeV. According to 
our model, however, one extra Higgs SU{2)l doublet, 
namely {Ha)\''^ , should have a mass much below M/, as 
was discusses in section |H] The small VEV found for it, 
in section IThI supported this view. Let us thus take its 



mass v{H 



A2 , 



Ma = Ml 



4 • 10* TeV, which is far 



above the lowest allowed value (~ 500 TeV) and does not 
lead to flavor changing neutral currents. With this value, 
the second term in l|VII.6|l leads only to a slight increase 
of Mj: M/ ~ 1.5 • 10^^ GeV. In general, the value of Mj is 
rather stable with respect to modifications of our model 
concerning the Higgs sectors i/(3,3, 1) and 77^(3,3,1). 
It is highly interesting that the value obtained for Mj is 
close (sec 1^) to the phenomenologically obtained mass 
scale (AsF^2,2}^ necessary to describe the mass squared 
difference observed in atmospheric neutrino oscillations. 
Morever, the same scale also describes the breaking point 
of the left-right symmetry. 

For the precise calculation of aL(M/) = Q!/j(M/) = 
cti,2{Mi) from HVII.5p . we need input masses for the D 
quarks and the leptons L. The mass of the 3'''^ generation 
D quark we take is based on the discussion about an ap- 
proximate left-right symmetry in the H and Ha sectors. 
We use 



](G^^) 



33 



4-10'' TeV 



(VII.7) 



Before renormalization, the lepton Lj, has the same mass. 
The ratios for the generation splitting of these quarks 
and leptons are cr* : cr^ : 1. The corresponding input in 
eq. I|VII.5|I allows now to calculate the values 0:3 (M/) 
and Q!i(Af/) — a2{Mj), which can then be used as initial 
conditions to go up to Mqut- After the study of the 
Yukawa coupling unification at Mgut, one can go back 
to the scales of the D and L states to find renormalized 
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values for their masses (see the next section and the ap- 
pendix). The corresponding change of eq. (|V11.5() will 
little affect the values of as and ai = a2 at Mj, from 
which one can start again. The result is 

a^^Mi) = ac\Mi) ~ 31.43 , 
a^\Mi) = a^\Mi) = al]j^{Mi) = 35.63 (VII.8) 

The D and L masses, found this way, are quoted in sec- 
tion IVII and have already been used in form of the mass 
matrices Md^ ej^G^^ and Me = elG^^ in sections Iml 
andHvl 

Above the scale of Mi G333 is unbroken and the quark- 
lepton states are unified together with the D states in 
(5l(3,1,3), (5ii:(l,3, 3) and the leptons L in L(3, 3, 1) 
multiplets. For the fermion masses we needed besides 
the VEVs from if(3,3, 1) also those from 77^(3,3,1). 
We take the masses of these Higgses to be negligeable 
for scales above M/ [similar to the mass of -ff(3, 3, 1)]. In 
fact, we have to do that because some members lie below 
Mi and the full [SU (3))^ symmetry must hold above M/. 
The corresponding ^-factors for ^ - Mj are therefore 

{bL,bR, bcf" = (-4, -4, -5) . (VII.9) 

With these values the meeting point gh — gn — 9c would 
be above the Planck scale because 6^ = fe^j is not much 
different from be- We know, however, from our treatment 
of the charged lepton sector, that the vacuum expectation 
values of 7?^ (6, 3,1) and 7?^(3,6, 1) play an important 
role. Since lying above the M/ scale, the masses of these 
two Higgses are equal due to the left-right Vlr symme- 
try: M(6,3, 1) = M(3,6, 1) = Mg. They contribute to 
the renormalization with the 6-factors 

(6L,fofl„6c)'= Q,^,o) . (VII.IO) 
We now have for /i ^ M/ 

ac\f^) = a~\Mj) - ^In^ , (VII.U) 

and 

-^^(A*-M6)%ln^ . (VII. 12) 

The Grand Unification Energy Mqut can now be ob- 
tained by setting /i = Mgut and equating (|V11.11() and 
(|VIL12|I . AfcuT depends on Me and increases with in- 
creasing Afg. It is interesting, that even for low values 
of il/g close to Ml we get a large values for Mqut- For 
instance for Mg ~ 3M/ we have Mqut ^ 10^^ GeV. Al- 
ready for Me- 5- 10^^ GeV we get Mqut ^ 3 • lO^^ GeV. 
Therefore, in our model we have Mgut ~ lO""^^ GeV, 
which thus insures proton stability compatible with 
present experimental limits. But we still have to see 
which restrictions are forced on us by top-bottom-tau 
unification. 



B. Top-bottom-tau unification 

In this section we study the running of the Yukawa 
couplings and their unification. We concentrate on the 
unification of the third generation couplings At, Ab, A,- for 
the top, bottom and tau fermions, respectively. In the 
SM, because of the small mixings in the quark sector, 
their evolution is little affected by the other couplings. 
In the considered model, the situation is different. Apart 
from the fermion couplings to H(3, 3, 1) [first coupling in 
(|1.8|l ] also couplings with Ha are important. In partic- 
ular, the Higgses iJ^(6, 3,1), i?^(3, 6,1) with common 
mass Ajfg < Mqut are important for gauge coupling uni- 
fication. Therefore, above the scale Mj, the following 
Yukawa couplings are relevant for renormalization: 

QlGqQr H + \LGlL H + QlAQQr Ha 

+\LA^L H\ + \LA^L H\ . (VH.13) 

We have to distinguish the coupling matrices Gg, Gl, 
A'^, A^, but have A^ = A^ due to the left-right Vlr 
symmetry which holds above /i = Mi. The elements 
of the diagonal matrices Gq, Gl determine the masses 
Md., Ml- respectively. 

As a consequence of the first term in (|VII.13|) one has 
already at the G333 level top-bottom unification: At(/i) 
and Ab(/i) must unify at Mi and evolve then further as a 
single coupling (fi) . This coupling should then unify 
with Ar(Ai) = ^L3{^^) at /I = Mqut. 

Below the scale Mi the coupling matrices Gq, Gl, 
A^ and A'^ = A^ split into more matrices depending on 
the Higgs field components they are attached to. In an 
obvious notation we have 

Q ^ I (_r , (_r , I , 

A^ ^ (A""'"^ , A^'""* , A"'^^^ . (VH.14) 

We left out the matrices A^^ , A^ and additional 
matrices from the neutral lepton sector. They are mul- 
tiplied with VEVs which are -in our model- small com- 
pared to competing terms in the same channel. In the 
approximations we use for the renormalization the G ma- 
trices remain diagonal and the diagonal elements of the 
matrices A remain zero. Furthermore, the matrices con- 
nected to A^ are the same as the ones from A^. But 
the matrices derived from A^ = A'" are no more strictly 
antisymmetric. 

The most important elements of the matrices (|VII.14p 
are the (3, 3) elements of the G's and the (2, 3) and (3, 2) 
elements of the A's: (G™)33 = At(^), (G'="'="')33 = A^(/i) 
etc. For the matrix elements of A'''^ we define 

{A''^^ = iAA(/i) , {A''%2 = -iAA(M) ■ (VH.15) 
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Clearly, we have Aa(m) = ^a(m) for fi - Mi. 

There is a restriction from the mass of the vector bo- 
son W for a combination of the VEVs multiplying the 
coupling matrices. With the notation 



e\ — vo sin (3 



— Vo cos f3 , 



(VII. 16) 



the condition is 



+ (/!)' + (/I)' + (/{l,3}3)^ 



(/ 



2(13} 



/{13}2 



(174 GeV)2 . (VII. 17) 



Since e\ and e| contribute to the masses of the third 
generation, the Vq term should be the dominant one. At 
the scale /i — Mz one has 



MMz) 



Vq sin /3 ' 

K{Mz) 



XbiMz 



Do COS 13 



Vq cos f3 



Here m° and m° are a little smaller than mf, and rrir 
respectively, since they refer to the diagonal parts of the 
down quark and charged lepton mass matrices. In sec- 
tions Iffl] and ^3 we found m^/rrib ~ 0.989, m^/rrir — 
0.966. 



We can now set up the renormalization group equa- 
tions for At, Ab, At, a a and A^- They are connected 
with each other and - due to the SU{?>)l x SU{?>)!i sym- 
metry dX fjL^ Mi - no other coupling intervenes. Below 

1^2 ^2 ^2 /I \2 

Ml we have ior -qt — , Vb ~ — ~ — " ■ — 



4tt 



Vr 



47r 



(VII. 18) and^^ 



(A/ 



VA 



4-R 



^ / 9 2 1 



17 9 \ 1 

—ai + -a2 -f 8a3 j -I- 9{^i - MA)-VtflA , 



(VII. 19) 



27r?7b = ^Vb + ^Vbilt + VbVr ~ Vb Q"! + ^"2 + 803^ + e{fj. - MA)vb + Va^ , 



(VII.20) 



(VII.21) 



1/1 \ /9 3 3 

27I"?7a = i^VA ivb + Vt- gtti - 3q;2 - 16q;3 j + 0{fi - Ma)va ( -^^Va + if] a ~ ^"i - ^"2 



(VII.22) 



1 3 A /9 3 3 

27r77^' = ( - — ai - -0:2 - 8^3 j + 6I(^ - MA)ry^ f -^^ + 3*7^ - - ^"2 



r 



(VII.23) 



At /i = Af/ the matching 

7?t = 776 EE 77Q3 , -qr = ?7L3 , f]A^flA = VA, (VII.24) 



2^'7a = 977^ + -VAVQ, - Va8 {aL,R + ac) , (VII.27) 



is required. 



Above Mi we have for rjQ^ , 77^3 , 77^ and 77AL - 4^ 
the equations 

2'^'7q, = 677^3 + r]Q^r]L3 + ^VQsVA - VQs^i^L.R + ac) , 

(VII.25) 

27r77^3 = 2r]l^ + 3r]L:,riQ^ - VL^^aL.R 

+e{fi - M6)3ru,VAL , (VII.26) 



2^'7al = VAL {}VL3 - 16aL,_R,) + 

0(m - M6)77AL (477AL - f aL.fl) . (VII.28) 

iA^ is the (2, 3) element of — and is only needed 
above M/. The matching condition at Mgut for the final 
unification of the couplings reads 

VQsi^IcvT) = VLsiMovr) , 7;^(AfGUT) = Val{Mgvt) ■ 

(VII.29) 
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'Concorde' 




2 4 6 8 10 12 14 16 
Logio[/i/GeV] 



FIG. 1: Unification of gauge couplings. Mi ~ 1.5-10" GeV, 
Me ~ 1.6 ■ lO^** GeV, Mgut ^ 2 • lO" GeV and a^^ ~ 39. 

The procedure of finding a solution with gauge and 
top-bottom-tau unification is the following: A given value 
of Mqut ~ 10^^ GeV (otherwise no solution is possible) 
fixes Mq. Taking then trial values for riQ^{MQ\jT) and 
77a(Mgut) and solving eqs. IIVII.25|) - (|VIlT28|l gives their 
values at Mj. These values determine rjt{Mj) = j]i,{Mj), 
rfriMi) and r]j^{Mi) = f}A{Mi). The rcnornialization 
group equations IIVII.19|) - (|VII.23p allow then to calcu- 
late Xt{Mz), Xb{Mz), Xt{Mz)- Clearly, the input values 
77Q3 (AfcuT) and riA{MQu^) have now to be changed such 
that Xr/Xb becomes equal to m^/m^ and At, Xb are in 
the perturbative region i.e. < 3. If this can be achieved, 
one can calculate from HV11.18|I Uq and tan (3 

«o = Tr + 4l^, tan/3 = (VII.30) 

Of course, only solutions with vq < v ^ 174 GeV are 
acceptable. 

C. Numerical solution for 

Mgut = 2 ■ lO'^ GeV, Md^ = Ma ^ 'i ■ 10^ TeV 

Here we present a numerical solution of the problem 
of gauge and Yukawa coupling unification in Eq, which 
satisfies all above mentioned requirements. We choose 
the unification scale to be 2 • 10^^ GeV, the masses of the 
heaviest D state and the Higgs field (77^)3'^ both equal 
to 3.98 • 10'' GeV. 

Further imput values are the third generation masses 

mt{Mz) = 173 GeV , toO(Mz) = 2.859 GeV , 

m°{Mz) = 1.689 GeV , (VII.31) 

the three gauge coupling constants ai fj, = Mz and a 
suitable value for r]t.b,T at A/qut 

Vt At (Mgut) ^ 0.0381 . (VII.32) 



'Bermuda Triangle' 



1 1 1 1 1 


1 1 l_ 






1 1 1 1 1 


1 1 1 


2 4 6 8 10 12 
Logio[/i/GeV] 


14 16 18 


: t ~ b ~ T unification: Ai(AfGUT) 


= XtiMavT) -- 



A,(Mgut) ^ 0.692. 

For this latter value all couplings remain in the pertur- 
bative region and vq < 1; = 174 GeV. As a result we find 
the solution 

Mg ~ 1.6 • IOI6 GeV , Mi ~ 1.5 • lO^^ GeV , 

?m(Mgut) = val{Mgut) ^ 0.0336 , (VII.33) 
with the following consequences: 

(Mgut) =^ 38.99 , 
r?t,fc(M/)=77Q3(M/)~ 0.0412, 
Vr{Mi) = 7]LAMi)^0M36 , 
VAiMi) = r]M) = Va{Mi) ~ 0.0368 , 

r?AL(M/) ~ 0.0598 , (VII.34) 

Xr{Mz) = Afc(Afz)^ ^ 0.612 , Xt{Mz) ^ 1.127 , 
vo ~ 153.48 GeV , tan^ ~ 55.59 . (VII.35) 

From the value found for wg, we can now determine 
(/I)' + (/{13}3)' from mrm. using then itTTTTm. 
(ITV^ together with el{G°^)33 = 3.98 • 10^ GcV and 
l/ll = \fil l/{i3}3l = we finally get 

/I = ±43.504 GeV , /| = T43.504 GeV , 

/{i3}3 = /^^"^ = 69.484 GeV . (VII.36) 

The solution for the gauge coupling and Yukawa coupling 
unification given here has been applied in the previous 
sections, in particular, for the evaluation of the renormal- 
ization parameters for all the different mass matrices. 

In figure 1 -'Concorde'- we show the evolution of 
the gauge couplings and their unification. Figure 2 - 
'Bermuda triangle'- exhibits the running of the Yukawa 
couplings r]t, Tib, Vt and their unification. In figure 3 - 
'desert spider' - the running of the (2,3) and (3,2) el- 
ements of the A-matrices and their unification is pre- 
sented. In these evaluations the splittings between the 
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'Desert Spider' 




Logio[/x/GeV] 



FIG. 3; Unification of (2,3), (3,2) elements of A*^ and 
matrices. .4J3(Mgut) = ^23(A^gut) 0.65 i. 



masses Mjj^ and M^. (which we discuss in an appendix) 
have been taken into account. 



it determines the mass scales for the light and heavy neu- 
trinos. We achieved a solution of the gauge and Yukawa 
coupling unification with strongly constraint parameters. 
It describes the evolution and the final convergence of 
many coupling matrices which differ significantly at low 
energies. The solution allows to calculate quite a num- 
ber of properties such as transition matrices from heavy 
to light fermions, Majorana phases and the double (3- 
decay matrix element. Due to the high unification scale 
(> 10^^ GeV), the model adequately suppresses dimen- 
sion six operators which induce nucleon decays. The pro- 
ton lifetime is above the presently accessible range. 

The presented iSg model can be supersymmetrized 
without changing the construction of the Yukawa sec- 
tor. A supersymmetric version would, however, affect 
the coupling unification picture given here. 
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VIII. CONCLUSIONS 

The Eq model presented has many attractive features. 
Only few input data are sufficient to obtain a realistic 
picture of the fermion masses and their mixings. The 
presence of new heavy fermions in the 'desert' plays an 
important role even for the mass matrices of the SM par- 
ticles. All generation mixings and CV violations arise 
from a single antisymmetric matrix A, which mixes the 
light fermions but also the light with the heavy fermions. 
The latter effect also contributes in an important way to 
the eigenvalues of the quark and lepton mass matrices. 
For instance, the main part of the /i meson mass and 
of the strange quark mass is generated by virtual transi- 
tions to heavy fermions. As a side remark we note, that 
the antisymmetric generation mixing matrix found here 
could lead to significant effects in rare weak decay pro- 
cesses with fixed phases of the new contributions. The 
matrix A, in combination with G, is also responsible for 
the bilarge mixing of the light neutrinos and their oscil- 
lation pattern. In the limit of no renormalization effects, 
the neutrino mixing is bimaximal. 

Those heavy new particles, which form 10-plets with 
respect to 50(10), have a hierarchical spectrum similar 
to the spectrum of the up quarks. The lightest ones are 
expected to lie in the low TeV region. 

The group Eq provides new insights about the unifi- 
cation of the three gauge couplings and about the unifi- 
cation of the Yukawa couplings of top, bottom and tau. 
The intermediate symmetry SU{'i)L x SU{i)R x S'C/(3)c 
with a discrete left-right symmetry plays a decisive role. 
The breaking point of this intermediate symmetry is fixed 
by the known gauge couplings gi and 52- Simultaneously 



APPENDIX: Renormalization Analysis 

Al. The D and L masses 

In our model the D and L states play a crucial role 
for the light families. Also, the mass splitting between 
these states must be taken into account in the study of 
the gauge coupling unification. So let's start with the 
renormalization of their masses. The Yukawa interac- 
tions occuring in ljVII.13|) obey the boundary condition 
G = Gq = Gl at ^ — Mgut due to Eq symmetry. We 
write 



Gq = Diag (AqjCT"* , - Aq^o-^ , 



and 



Gl = Diag(ALiCr'* , - Al^ct^ , Xl,) 



(A.l) 



(A.2) 



For a given value for Ag — Al at A/qut the values for 
Aq3, Al 3 at ^ = Mj can be obtained using eqs. HVII.25|) - 
(|VIL28p . The gauge interactions contribute to the run- 
ning of Aqjj2) s-nd Al^2) ^ similar way as for Xq.^ and 
Al;,. However, also the Yukawa interactions are to be 
taken into account. One has 



Mi 



^Q2 
^Q3 



Ac 



Mi 



A 



pi 



Mgl 



Mi 



Aio 



Mr 



Mgut 



(A.3) 



(A.4) 
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where the P-factors are defined as follows 



Pal = cxp 



(A.5) 



All numerical analysis are carried out according to the 
solution presented in section IVll CI The P-factors have 
the values 



Pq, ~ 1.030 , Pl, ^ 1.034 , 
Pa 1.027 , Pal ^ 1-015 . (A.6) 

At the scale Mj the matrices Gq and Gl are given by 
Gq{Mi) = Diag (^Kfa\ - ^a^, l) Xq^Mi) , (A.7) 



i_4 



Gl(M/) = Diag (<t7 
where 



4^2, 1)Xl,{Mi) , (A., 



1 p3 p3 



3 



^2 = 



(A.9) 



The coefficients k„ and the factors of cr'* and 

occuring in G at the GUT scale [see (IA.31II ] as obtained 
from mu at ^ = Mz [eq. (|III.4p ]. (|A.7p and (jA.Sp allow 
to determine at Mj the ratios M^./Md^ and MlJMl^ 
since M^i ~ Gq, Ml Gl- Below M/ the corresponding 
ratios run due to gauge interactions. We have 



Mb 



Mr 



A 
X 



Mn 



Mr 



X pl{' {Md,)pI, {Md,)p'o.I" [Md, )p^l [Md, ) (A. 10) 



The scaling factors pi are defined in ljA.25|) in a similar 
way as the P-factors. In our model we find 

P^,{Md.,) ^ 1.034 , Pc.,{Md^) ^ 1.066 , 
p„,(Afz33) ~ 1.025 , p„3(Md3) ~ 1.040 . (A.ll) 

The ratios (IA.10|) do not run below the scales p = Mdi 
and Md^ respectively. Combining (|A.7p and ljA.10|) . we 

get for the mass ratios of the Z?-statcs 

M^^ = e^G^^ = Diag - K^^a\ l) Md, , 

(A.12) 



1.237 , 



0.965 



(A.13) 



For the mass ratios of heavy quarks to heavy leptons 
we have a,t p = Mi 



Ml 



Mi 

^Q3 



p3 p3 p— 1 p— 3 
^Qs^A^L, ^AL ' 



Mr 



Mr 



Pn P7^ 

Q3 L3 



(A.14) 



Using HVII.34|) and HA.5p . (|A.14|) gives 



M, 



Di 



Mr 



= 0.960 



Mr 



M 



L2 



0.928 . {A.15) 



Mr 



These values form the starting point for the further run- 
ning of these mass ratios down to their own mass scale 
due to their gauge interactions 



Mr 



Ml^ 



Md, 



Ml. 



n 



(A.16) 



ba are the 5-factors in the interval Md^ — Mj and 

V5' ' / ' a VlO'2'"/ ' 

a=iUil)Y,SUi2)L,SU{3)c) ■ 



(A.17) 



The heavy quark to heavy lepton mass ratios at their 
mass scales [see (|VII.7ll . (|A.12|I . (|A.13ll ] are found to be 



^^^^ - 1 572 ^'^^^ 



1.249 



Mr 



M 



L3 



1.039. (A.18) 



These mass splittings do only minimally affect the gauge 
coupling unification. Their effect is a subleading one only. 
By combining ljA.13|) and l|A.18p . one can find the mass 
ratios for the P-states 



M 



LL 



with 



^G^^ = Diag (/^f^^7^ - l) Ml, , 

(A.19) 



Kf ^ ~ 0.817 



K2 



0.802. 



(A.20) 



Recall that in the application of our model we use Md, = 
3.98 • 10^ GeV. 



A2. The running of the matrices G"", 
G''^ and G""^^ 

Below the scale Mj instead of the matrix Gq we have 

two matrices: G"" and G'^'^. G"" accounts for the up 
quark Yukawa couplings. At p — Mz it is given by 
(|III.3p . (|III.4|) . Its hierarchical structure is dictated from 
the values of up- type quark masses pi.lO|l . The diago- 
nal matrix G'^'^ describes the coupling of qd with H^'^. 
The matrix Gl is also splitted below p = Mj, namely 
to G*^ and G"'', which describe the couplings le^H^^^ 
and WHI'^, respectively. At the scale Mj these couplings 
satisfy the conditions 

G""(M,) = G^^Mi) = Gq{Mi) , 



G" 



" {Mi) = G''<'{Mi) = Gl{Mi) . (A.21) 



Since we have explicit information on G"" at Mz, we can 
calculate G at the GUT scale and then derive G'^'^, G^ 
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at /X = Mz- G^^ is connected with the neutrino sector, 
the renormahzation of which is studied in appendix A4. 

Below Ml the running of the up and charm quark 
Yukawa coupHngs are due to foUowing equations 



17 9 

47rA^ = A„ ( 3?7t - — ai - -a2 - 8013 



(A.22) 



/ 17 9 1_ \ 

47rA'^ = Ac ( 377t - — ai - -a2 - Sas + 6l(/i - Ma)^Va j > 

(A.23) 

It is convenient to introduce for each coupUng ?/i 
(?7t, Qfi, q;2, • • • ) the quantity 



Pi(p) = exp 



1 



47r 



(A.24) 



These scahng factors aUow to express the couphng con- 
stants at arbitrary scales and thus are useful for satisfy- 
ing the matching conditions on boundaries. According 
to IIA.241). we also have 



exp 



1 

47r 



rii{p')d\n pt' 



= P»(mi)M(M2) . (A.25) 



Taking all this into account, equations (|VII.19|I . (|A.22|I . 
(|X23)l give 



Au 

a7 



Mz 



At 



p'/'{Mz)pl/'{Mz)pli'{MA) , (A.26) 



Our model gives 



Ku 1.528 , Kc ~ 1.386 



(A.32) 



From ljA.31|) , it is easy to derive G''''^ and G'^ at the 
scale Mz- The result is 

G^^Mz) - Diag («f>, - Kfa\ l) • \,{Mz) , 

(A.33) 

G^'^\Mz) = Diag (nl^W - kI^W l)-\r(Mz) , 

(A.34) 

^ p,{Mz)Pt\Mz)p^,{MA) , 
Atf = p,(A/z)pr'(Mz)p^jMA)p^i(MA) , (A.35) 



with 



5+ _ 1 



(A.36) 



The numerical values of these factors are 



Ae 

At 



Mz 



At 



p'/\Mz)p',''iMz)p;,':{MA)p=^''iMA) 



1/2 



1/2/ 



-1/2. 



(A.27) 

At the scale Mi the ratios \u,c/M can be expressed by 



^Q3 

All 

a7 



a/gut 



according to (jA.3|l . Therefore, we have 



Au 



Ma 



At 



Mz 



Xc 

At 



Ma 



(A.28) 



where 

K„ = P^^Plp'J'{Mz)pl/'iMz)pl['iMA) , 

(A.29) 

Numerically wc find in our specific model 

Pt(M2) ~ 1.130 , p6(M2) ~ 1.112 , 

Pf,jMA) ^ PnA^A) ^ 1.040 . (A.30) 

Consequently, if G has at p — Mz the form given in 
(jIII.4|l . we have to take at the GUT scale 



C?(Mgut) = Diag ( —a\ - — ct^, 1 ) • At(M, 



GUTj ■ 

(A.31) 



Kf = 1.023 , 
f e+ = 0.804 , 



"-2 



0.984 , 
= 0.848 . 



(A.37) 



A3. Renormalization of , A'" matrix elements 



At the scale A/qut the antisymmetric matrix A, which 
describes the fermion couplings with the components of 
the Higgs field Ha, is postulated to have the form 



A(Af, 



GUTj 



, 


icr , 


— icr 


-icr , 


, 


75 


icT , 


"73 ' 






XaV2. (A-38) 



Between the scales Mi and Mqut , instead of one matrix 
A we are dealing with the two matrices A^ and A^ = A^ 
in the Yukawa interaction of l|VII.13|) . Due to RG effects, 
they will differ from the original matrix A{Mq\j'j^). A'^ 
and A^ remain antisymmetric above Af/ and have the 
forms: 
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, 



, 



1 

\/2 



, 
]_ 

V2 



1 





(A.39) 



Here k*^'-^ and Aa,A^ are scale dependent renormaliza- 
tion factors, which can be determined through the RG 
equations of the A'^ and A'^ matrices: 

167r2(AO)' = SAQAQ^AQ + IAqA^qAQ + ^AQXqX^q 
+AQ{3TTiAQAQ^)-8gl^-8gl) , (A.40) 



A^'^iMz) 




-1.2a 

(A.46) 



e{n^Ah)A^ {iA^^ A^+\Tr{A^ A^^y^gl j^) (A.41) 
At the scale Mqut, the boundary conditions are 



GUT J 



A^{KU 



GUT J 



and 



(A.42) 

Fro m lfO0)l and (fOT)) follow the RG equations HVII.27|) 
and l|VII.28|) for Xa and A^ , which we solved numerically. 

It is easy to observe, that the ratios A^g/A^g and ^13/^^3 
do not run, while the other ratios are scale dependent. 
For the factors k*^^^ we have at /i = Mj 



,3/2 



1/2 



(A.43) 



where Pq^, Pl3 are given in HA.5|I . (|A.6p . 

Below the scale Mj, instead of one A''^ matrix there 
are matrices which represent the couplings of colored 
fermions with Higgs doublets and singlet. Namely: 



Ha ^ A'^'^qd {Ha)2 



A'^^qD (Ha) 



1.2 



A^'^Dd {Ha)1 

(A.44) 



At the scale M/ all matrices unify ^''''(Af/) = 
A°^{Mi) = A'^^{Mi) = A^{Mi). There they are pre- 
cisely antisymmetric and their matrix elements match 
with the appropriate k^j, factors of (|A.39p . RG study al- 
low to calculate these matrices at the scale needed. Our 
numerical analysis gives the following results for the A'^- 
matrices at the scale = Mz (where contact with exper- 
imental data can be performed): 



A'^'^{Mz) 



, 
-1.16cr 

0.99cr , 



1.16ct , — cr 

0, 75 \Xa{Mz)V2, 




0.99 



(A.45) 



The couplings are 

Xa{Mz) = 1.262 




X\{Mz)V2. 



(A.47) 



Xa{Mz) = 1.131 



A5i(Mz) = 1.220 . 



(A.48) 

A similar analysis can be performed to obtain the gen- 
eration mixing matrices A^ for leptons. Below Mi we 
are dealing with three types of matrices , A^ ^ 

and A^ ^ . These coupling matrices are relevant for the 
charged lepton sector. At the scale n — Mz we find 



A"" 



'{Mz) 



, l.llcr , -a \ 
-1-11^, 0, ^ A^(A^^)^/2, 
1.04a , , y 

(A.49) 



A' 



'{Mz) 



, 1.14(7 , 
-l.llcr , , 
1.06a, -J^, 




A" 



'{Mz) 




The coupling factors are 



X\{Mz)V2, 



(A.50) 



Xf{Mz)V2. 



(A.51) 



Aa(A/z) 



= 0.957 , 
A!^(Mz) 



AS(Mz) 
= 1.068 . 



0.923 



(A.52) 
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The picture which shows the running of all (2, 3) and 
(3, 2) elements of the matrices A'^ and and their final 
unification is presented in figure 3, the 'desert spider'. 



A4. Neutrino mass matrix renormalization 

Here we will perform the renormalization analysis for 
the neutrino sector. 

Dimension five operators, responsible for neutrino 
masses, are generated by integrating out the 'right 



handed' states (Lf) 



2Ja 



The 



masses are deter- 



mined by the matrix F^^'^^S. The matrix S describes 
the generation dependent Majorana couplings of these 
states to the symmetric sextet component of the Higgs 
field Hs- In our model S is postulated to be the bihnear 
matrix product in generation space (jV.Sp . We take this 
form to be valid at Mi with G ^ Gl and A^ A^. With 
the appropriate scaling factors, at fi = Mi the matrix S 
has the form: 



1 



^2 K^ct'^X 



S/Xs 



4 

-ax 



(A.53) 

Again, as in ljV.9ll , only the leading terms in a are exhib- 
ited here. The renormalization coefficients k are shown 
in (|A.75|) . At the scale M/, the mass matrix for the states 
i>„ is 



a/? 



(A.54) 



The eigenvalues are fii ~ 

fi2 ~ a^i4f\sF^^'^^ and \sF^^'^\ As we can 
see, these scales are separated by large distances. 
Because of this fact strong renormalization effects occur. 
They are the cause of the difference between the neutrino 
mass matrices HV.10|) and (|V.13ll . The decoupling of 
the three i)a states occurs step by step. Thus the 
renormalization has to be performed separately in each 
ener gy i nterval p^ . [l5l |. By generalizing the results of 
ref. flql we present model independent formula for the 
running of the neutrino mass matrix mi, and then apply 
them to our model. 

The couplings in the neutrino sector involve Dirac and 
Majorana mass terms. One can choose a basis in which 
the mass matrix for the heavy neutrinos is diagonal. 
Thus, without loss of generality, one can write the cou- 
pling terms in the form 



(A.55) 



The matrix A '^'^ is related to the original Dirac matrix 

^0 



A 



Dirac 



\ Dirac t jT 



(A.56) 



Here the unitary matrix Us diagonalizes the matrix M^o: 



(A.57) 



By integrating out the state Va, the light neutrino mass 
matrix gets a contribution at the scale /Zq. Without 
renormalization effects, rrii, would have the form 



m 



a/3 



-{Yi+Y2^Y:irP 



(A.58) 



where 



Y°'P _ J_^^Dirac^Q2^^Dirac'j/32 

Y;^'^ = J_(^Dirac)a3(^Dirac)/33 (^.59) 



The division of the mass matrix in three parts is conve- 
nient in order to see the contributions coming from each 
integrated state Va- Each d = 5 operator (Ki), gener- 
ated on the scale /i^, runs from this scale down to Mz 
according to the RG equations 

4^|i^"^ = i;"''(677t - 3a2) , with (a, /3) ^ (3, 3) , 

(A.60) 



4^^r,3" = Yi"{Qi^t + - 3a2) , with a ^ 3 , 



47r-y/3^r,33(6r;t + r;,-3a2) ■ 
These equations have the solutions 

Yt^{^^')^Yt\^i)r,{n')r-\^^) , 



(A.61) 



(A.62) 



(A.63) 



y/"(M') = Y?-{^i)r,{^,')rr{^i')r-\^^)r-\^i) . (A.64) 



Y,^H^^') = Y,^^^i)r,i^^'yMr;\^,)r-'{^,) , (A.65) 
where 

rgi^i) - pi'{fi)pl^{n) , rr{^^) = p-'/'ifi) . {AM) 

Before the emergence of the d = 5 operators, only the 
Dirac couplings run. They obey the equations .24 | 



Dirac \ a/3 



= 



(a,/3)^(3,3) , (A.67) 



(A.68) 
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(A.70) 



According to these equations, we have for /i < Mj trino mass matrix at the scale /i = Mz 

(;^Dhac)3a(^) = (AD"-^)^" (Af,)r, (/i) . (A.69) " ^ ^ ' ' 

With all these results, one can write down the light neu- where 
I 



5^H^^rl{Mz) [Yi^ + r-Hp2)Yi^ + r-\ni)Y^^) . (A.71) 



The quantities Yi are given in l|A.59|l and the renormal- 
ization factors in HA.66|I . 

We now apply this result to our model. The matrix 
j^^Dirac-j built according to HA. 56(1 taking 

A?"'^^(A/,) = Diag {n[a\ - 4^^^ l) A,(A//)e}(M,) . 

(A.72) 

The value of e\ at the scale Mi can be calculated from 
the RG equation 

4^^e} = e\ (^-3?7t + + -^a^ , (A.73) 

knowing its value at the scale /i = Mz'- 

e\{Mi) = e\{Mz)pi\Mz)p'J^\Mz)p%\Mz) . 

(A.74) 

In our model, it turns out that eJ(M/) ~ 124.5 2 GeV. 
The values of the factors appearing in ljA.53|) are 

~ 0.707 , K2 0.746 , ~ 1.017 . (A.75) 



After constructing the matrix A^"''*'^, one can calculate all 
elements of Yi according to ljA.59|l . The numerical val- 
ues of the renormalization factors, appearing in (|A.70p . 
I|A.71II are 

rg(Mz)~0.58, r<,(A*i) -0.75 , 

7-3(^2) ^ 0.87 , r,{Mz) ~ 0.96 . (A.76) 

The diagonalization of the neutrino mass matrix ob- 
tained in this way allows to calculate the neutrino mixing 
and the ratio R = Am'^^i/ Am'^^^ from the single param- 
eter X. To obtain values for R and the mixing angles, 
which lie within experimental bounds, we have to use 
X ~ 2.8 which is bit smaller than the value x — 3.5 found 
without renormalization corrections !25l |. The mixing is 
now no more bimaximal, but still bilarge. The results 
are discussed in the text (section O. 
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